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Abstract 

We consider photon radiation in it events at the upgraded Fermilab Tevatron 
and the CERN Large Hadron Collider (LHC) as a tool to measure the electric 
charge of the top quark. We analyze the contributions of itj production and 
radiative top quark decays to pp — > "yi^ubbjj, assuming that both 6-quarks 
are tagged. With 20 fb _1 at the Tevatron, the possibility that the "top quark" 
discovered in Run I is actually an exotic charge —4/3 quark can be ruled out 
at the ~ 95% confidence level. At the LHC, it will be possible to determine 
the charge of the top quark with an accuracy of about 10%. 
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I. INTRODUCTION 



It is widely believed that the new particle discovered at Fermilab in 1995 |!]|| is the 
long sought top quark. Its behavior is certainly consistent with Standard Model (SM) 
expectations for top, but many of its properties are still only poorly known. In particular, 
the top quark's electric charge — one of the most fundamental quantities characterizing a 
particle — has not been measured yet. It still remains not only to confirm that the discovered 
quark has charge 2/3 (and hence the expected SM quantum numbers), but also to measure 
the strength of its electromagnetic (EM) coupling to rule out anomalous contributions to its 
EM interactions. 

Indeed, one alternative interpretation has not yet been ruled out: that the new particle is 
a charge —4/3 quark. In the published top quark analyses of the CDF and D0 Collaborations 
(see Ref. || for a review), the correlations of the 6-quarks and the W bosons in pp — > it — > 
W + W~bb are not determined. As a result, there is a twofold ambiguity in the pairing of W 
bosons and 6-quarks, and, consequently, in the electric charge assignment of the "top quark" . 
Besides the Standard Model (SM) assignment, t — > W + b, "t" — > W~b is also conceivable, in 
which case the "top quark" would actually be an exotic quark with charge q = —4/3. 

Interpreting the particle found at Fermilab as a charge —4/3 quark is consistent with 
current precision electroweak data. Current Z — > £ + £~ and Z — > bb data can be fit with a 
top quark of mass m t pa 270 GeV, provided that the right-handed 6-quark mixes with the 
isospin +1/2 component of an exotic doublet of charge —1/3 and —4/3 quarks, (Qi, Q^r @- 
If the top quark would have a mass of m t ~ 270 GeV, it would have escaped detection in 
Run I of the Tevatron. In this scenario, the particle discovered in Run I is the Q4. A direct 
measurement of the top quark charge would provide a definitive test of this interpretation. 

There are several techniques to determine the electric charge of the top quark, q top , in 
future collider experiments. At a linear e + e~ collider, the full set of form factors describing 
the most general ytt vertex function compatible with Lorentz invariance can be probed 
in e + e~ — > it || (similarly for a muon collider). A photon-photon collider would provide 
even more precision via the two electromagnetic vertices in 77 — > it. The status of future 
lepton and photon colliders is unfortunately still somewhat uncertain, but in the meantime 
top quarks will be produced copiously at the Fermilab Tevatron and the CERN LHC. At 
these hadron colliders, it production is so dominated by the QCD processes qq — > g* — > it 
and gg — > it that a measurement of the jit form factors via qq — > 7* — > it is hopeless. 
Instead, one can measure the top charge by measuring the charges of its decay products, 
the final state 6-jets and W bosons. One can further attempt to determine not only the top 
quark's charge but also the EM coupling strength by directly studying top's electromagnetic 
interactions through photon radiation in it events. 

The first method — measuring the top charge by reconstructing the charges of its decay 
products - - is difficult to realize for it — > all jet decays; however for di-lepton decays, 
it — > £ + v£~vbb, and for semileptonic decays, it — > £ vjjbb, it should be feasible. The 
6-jet charge can be determined from a measurement of the charges associated with the 
tracks in the jet. A preliminary measurement of the 6-jet charge using Run I CDF data || 
shows a slight preference for the SM top charge assignment. The direct measurement of 
the 6-jet charge has the disadvantage that many tagged 6-quarks are needed to obtain a 
statistically significant result. In a given event, missing or out-of-cone tracks can distort the 
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measurement. Information on the 6-jet charge can also be obtained from the charge-sign 
of the lepton, £ b = e, //, in events where the 6-quark is identified via a soft lepton tag, i.e. 
where the 5-quark decays semileptonically, b — > ct^v. In the absence of B d — B d and B s — B s 
mixing, the charge-sign of the lepton is directly correlated with the charge-sign of the parent 
6-quark. The difficulties associated with using soft lepton tagged 6-quarks are low efficiency 
due to the small b — > civ branching ratio and the presence of wrong sign leptons originating 
from Bd — Bd and B s — B s mixing, which one has to correct for. In addition, if only the 
charge of the b decay lepton is measured, only the sign but not the magnitude of the charge 
of the 6-quark is determined. 

In this paper, we explore the possibility of measuring the electric charge of the top quark 
at hadron colliders through photon radiation in it events, assuming that the EM coupling 
strength takes its SM valueQ. Because top quarks can radiate photons in both top quark 
production and top decay, we consider both processes: 

pp — > it'y and pp —>it, t — > Why. (1) 

We concentrate on the lepton+jets modeQ, 

pp -> ^ubbjj, (2) 

£ = e, fi, and assume that both 6-quarks are tagged. We present results for the SM top 
quark and for comparison we also consider a charge —4/3 quark of the type discussed in [|J. 

Photon radiation in top quark events was not observed by the Tevatron experiments in 
Run I. In Run II, an integrated luminosity of 2— 20 fb _1 is envisioned, and a sufficient number 
of j^ubbjj events may be available to measure the top quark charge. At the LHC, where 
the cross section for it production is more than a factor 100 larger than at the Tevatron, 
one can hope for a determination of q top with a precision of about 10%. 

Details of our calculation are given in Sec. II. In Sec. Ill we present results of our numeri- 
cal studies. In addition to total cross sections, distributions for the transverse momentum of 
the photon and the separation between the photon and the 6-quarks in the pseudorapidity - 
azimuthal angle plane are found to be sensitive to q top . In Sec. IV we perform a quantitative 
analysis of how well the charge of the top quark can be measured at the Tevatron and the 
LHC. Finally, in Sec. V, we present our conclusions. 

II. DETAILS OF THE CALCULATION 

Our calculation is carried out at the tree level. We assume that both 6-quarks are tagged 
with a combined efficiency of 40%, which is included in all our numerical results. Top quark 



1 This is required by electromagnetic gauge invariance for on-shell photons radiated in tt events. 

2 Although di-lepton events, p'p — ► ^yl^ul^i/bb, are cleaner than lepton+jets events, their branch- 
ing fraction is about a factor 6 smaller than that of the lepton+jets mode. The all-jets mode 
has a larger branching ration than lepton+jets, but is plagued by a large QCD background and, 
therefore, is also not considered here. 
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and W boson decays are treated in the narrow width approximation. In this approximation, 
there are three contributions to pp — > / y£ ± ubbjj: 

• Radiation in top production: pp — > ttj — > / y£ ± vbbjj. 

• Radiative top decay: pp — > — > / yW + bW~b — > / y£ ± vbbjj. Here either the t or t quark 
decays radiatively. We denote these decays generically by i — ► W67. 

• Radiative decay: pp — > — > W + bW~b — > , y£ ± ubbjj. In this case, one of the W 
bosons decays radiatively, — > £2/7 or — > jj'7. 

Both tt7 production and radiative top decays are sensitive to q top . On the other hand, 
radiative W decays are not, and should be suppressed by cuts. This can be achieved by 
requiring that 

m{jh) > 90 GeV and m T (£ T , 1> T ) > 90 GeV, (3) 

where j) T denotes the missing transverse momentum originating from the neutrino in the 
event which is not observed, and m(jjj) is the invariant mass of the jjj system. mr(^7;pV) 
is the ^7j$ T cluster transverse mass, which is given by 

m 2 T (£rJ T ) = (^PWl)+m 2 (£l)+^ 2 - (p T (£ 1 )+^ T ) 2 , (4) 

where Pt(^7) and m{£ r y) are the transverse momentum and the invariant mass of the £j 
system, respectively. For W — > £2/7, the cluster transverse mass sharply peaks at the W 
boson mass, Myy- In the narrow width approximation, the cuts listed in Eq. (Q) completely 
eliminate the contributions from radiative W decays to the r y£ ± i>bbjj final state. We thus 
ignore radiative W decays in the following. 
The matrix elements for 

99, qq^ *<7 l^vbbjj (5) 

are calculated using MADGRAPH and the HELAS library ||, including the spin correlations 
for the subsequent top decays. To compute the matrix elements for 

gg, qq -»• it, t -> iy + &7, t -> W~b (6) 

and its charge conjugate sibling, we use the crossed form of the W — *> £67 matrix elements 
of Ref. ||. We ignore spin correlations in the decay of the W boson in t — > W + b r y. For the 
i — > W~b decay, spin correlations are taken into account. All numerical results presented 
below are obtained summing over electron and muon final states. 

The top quark charge, q top , is treated as a free parameter in our calculation. The electric 
charge of the 6-quark, q b , is related to q top by 

qb = qtop - qw, (7) 

where qw = ±1 is the charge of the W boson. The sign of qw depends on whether the 
third component of the weak isospin of the "top quark" is T 3 = +1/2, as in the SM, or 
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T3 = — 1/2 as in the scenario of Ref. 0. We assume the Wtb coupling to be of V — A form 
and BR(t — > Wb) = 1. The fermions originating from W decays are assumed to be massless. 

We note for future reference that the differential cross sections for gg — > it'y and for 
qq, gg — > it with one top quark decaying via t — *> Why are invariant under the transformation 

qtop — > —qtop, qw —> —qw, qb — ► — <?&• (8) 

Since interference effects between the incoming light quarks and the final state top quarks in 
qq — ► it'y are very small (they contribute less than 1% to the cross section), the cross section 
for pp — > ryi^vbbjj is to a very good approximation invariant under the transformation (^j). 
This will be relevant when we estimate how well q top can be measured in future Tevatron 
and LHC experiments (see Sec. IV). 

For all our numerical simulations we have chosen the set of SM input parameters to 



be: m t = 175 GeV ggjn), m b = 5 GeV, M w = 80.3 GeV, sin 2 6 W = 0.23 and a(M, 
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1/128 flT2" |. We employ MRSR2 parton distribution functions |T3]| , selecting the value of the 
factorization scale to be /if = Vs, where s is the squared parton center of mass energy. 
Calculations are carried out for pp collisions at y/s = 2 TeV (Tevatron) and pp collisions at 
v^= 14 TeV (LHC). 

To simulate detector response, we impose the following common transverse momentum, 
rapidity and separation cuts at the Tevatron and LHC: 

Pr(b) > 15 GeV, |t;(6)| < 2, (9) 

p T (£)>20GeV, \r}(£)\<2.5, (10) 

Pr(j 1,2) > 20 GeV, \ V (j 1<2 )\ < 2.5, (11) 

i> T > 20 GeV, &R{i,j) > 0.4 for i ^ j, (12) 

where i,j = 7, £, b, b, ji, j2- ji and j'2 denote the two jets originating from the hadronically 
decaying W in the event. As mentioned above, we assume that both b quarks are tagged 
(40% combined efficiency) so that b and b are distinguishable from ji and In addition, 
we require that 

Pr{l) > 10 GeV, \r)(j) \ < 2.5 at the Tevatron and (13) 
p T (l) > 30 GeV, |r/( 7 )| < 2.5 at the LHC. (14) 

Here, r\ is the pseudo-rapidity, and 

AR(i,j) = [(A$(*,j)) 2 + (A V (i,j)f] 1/2 (15) 

is the separation between two particles % and j in the pseudorapidity - azimuthal angle 
plane. 

The transverse momentum cuts listed in Eqs. (|ToD - ( [i~2]) and in Eq. ([i4[) for the LHC 
are sufficient for running at low luminosity, C = 10 33 cm _2 s _1 . For operation at the design 
luminosity of £ = 10 34 cm -2 s _1 , the px cuts must be increased due to pile-up effects caused 
by the large number of interactions per beam crossing. In the following we only consider 
^t^vbbjj production at the LHC in the low luminosity environment. 

We do not consider any potential background processes here. The dominant background 
is expected to originate from W7+ jets production and should be manageable in much the 
same way as the W+ jets background for ti production. 
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III. PHENOMENOLOGICAL RESULTS 



In this section we present our results for total cross sections and various photon dis- 
tributions. We look at contributions from individual subprocesses (radiative production, 
radiative decays) and compare results for the SM top quark to those for an exotic charge 
—4/3 top (as in Ref. in order to determine which quantities are sensitive to the value of 
Qtop- 

The events passing the cuts listed in Eqs. (|3|) and fllOD - flH| ) can be split into three 
separate samples, each designed to enhance one of the three subprocesses: (i) radiation in 
top production ("££7 selection cuts"), (ii) radiative top decay with leptonically decaying W 
("i — > Why — > iuby selection cuts"), and (hi) radiative top decay with hadronically decaying 
W ("t — > Wbj — > jj&7 selection cuts"). 

The event sample for radiation in top production (tt'j selection cuts) is obtained by 
suppressing radiative top decay events. We do this by selecting events which satisfy 

m{b 1>2 jn) > 190 GeV and mrfai^rM > 190 GeV - ( 16 ) 

Here, b\, b 2 = b,b and b\ ^ b 2 . We note that by imposing the cuts on either 6-quark 
combination, we conservatively assume that the charges of the individual b quarks have not 
been determined. 

The total cross sections for q top = 2/3 and q top = —4/3, imposing the itj selection cuts 
of Eq. (|i~6|) and the cuts described in Sec. II, are listed in Table |a for the Tevatron and the 
LHC. At the Tevatron, the itj cross section is completely dominated by qq annihilation. As 
a result, photon radiation off the initial state quarks constitutes an irreducible background 
which limits the sensitivity of the cross section to q top . In contrast, at the LHC more than 
75% of the itj cross section originates from gluon fusion. Since gluons do not radiate photons, 
the LHC cross section scales approximately with qf op . Radiative top decays contribute about 
10 - 20% (1%) to the cross section for q top = 2/3 (q top = —4/3) in this phase space region. 
At both the Tevatron and LHC, increasing |g top | increases the total cross section with tt'y 
selection cuts. 

The photon transverse momentum distributions for the individual contributions to the 
r y£ ± i , bbjj differential cross section at the Tevatron (LHC) for it'j selection cuts are shown in 
Fig. |l] (Fig. |) for both q top = 2/3 and q top = —4/3, keeping the same vertical scale for both 
charges. The cross section from it production followed by t — > Wb'j — > iuby is about 20% 
larger that from pp — > it followed by t — > Why — > jjtry over most of the pt range considered. 
At the LHC, the Pt{i) distribution for the radiative top decay channels drops much faster 
than that for pp — > itj. 

A sample for radiative top decay with leptonically decaying W (t — > Wb-f — > iuby 
selection cuts) is obtained by requiring that 

7^(61,2^7; J* r ) < 190 GeV and m (b 2 ,dn) > 190 GeV ( 17 ) 

in addition to ( P5p being not satisfied. These cuts enhance the contribution of the process 
pp — > ti, with t — > Wb^f — > Ivky. Similarly, to get events with radiative top decay and 
hadronically decaying W (t — > Why — > jjby selection cuts) we impose 

m T (61,2^=7;^) > 190 GeV and 160 GeV < m(b 2 ,xjjj) < 190 GeV (18) 
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and require that (|T^) is not fulfilled. We thereby obtain an event sample in which the 
process pp —> ti, with t — > Wb'y —> jjb'y is enhanced. The cuts are imposed on either of the 
6-quark combinations, i.e. we assume that it has not been determined whether the b or the 
b originates from the radiative top quark decay. Since the finite experimental resolution of 
detectors significantly broadens the Breit-Wigner resonance for t — > Why — > jjby, we have 
chosen a relatively large window for m(&2,iJJ7) around mt. Requiring that events do not 
satisfy ([16]) significantly improves the efficiency of the t — > Wb'y selection cuts. 

The cross sections for pp — > ^l^vbbjj imposing t — > Wb'y — > ivby (Eq. (|17D ) and 
t — > Wby — > jjb'y selection cuts (Eq. fllTf)) are shown in Table 0b and |I|c for q top = 2/3 
and q top = —4/3. The total cross sections are smaller than the corresponding numbers 
for tt'y selection cuts. Moreover, the cross section for it production followed by radiative t 
decay decreases by a factor 3 to 4 when the top quark charge is varied from q top = 2/3 to 
Qtop = —4/3. This is due to interference effects between the diagrams where the photon is 
radiated from the t-quark line and the diagrams where the photon is either emitted from 
the final state W or the 6-quark. For the photon p? cut used, a significant fraction of 
the events in the SM at the Tevatron in the regions defined by the t — > Wb'y — > ivby 
and t — > Wby — > jjb'y selection cuts originates from tt'y production. For q top = -4/3, 
tt'y production is in fact the dominant source in these regions. Nonetheless, the net result 
is a decrease in the total cross section when we switch from q top = 2/3 to q top = —4/3. 
Recall that, in contrast, the cross section increased with the change in top charge for tt'y 
selection cuts. This suggests that the ratio of the cross section associated with radiation in 
top production to the radiative decay cross section would be very sensitive to the top quark 
charge. This issue is explored in the next section. 

Combining the cross sections for all three phase space regions, we expect about 60 
ryl^vbbjj events at the Tevatron for 20 fb _1 . The number of events expected at the LHC is 
about an order of magnitude larger: 550 'yl^vbbjj events for 10 fb _1 . 

Returning to the radiative decay sample, we show in Fig. |3| the individual contributions 
to the r y£ ± vbbjj differential cross section for t — > Wb'y — > Ivby selection cuts at the Tevatron 
for q tov = 2/3 and q top = —4/3 as a function of the photon pt- The individual contributions 
at the LHC are shown in Fig. |. We see that the contribution from tt'y production in the SM 
at the Tevatron exceeds that from pp — > it with t —>■ Wb'y — > ivby only for Pt(j) < 25 GeV 
(Fig. |3|a). Increasing the photon p? cut thus improves the efficiency of the t — > Why — > 
Ivby selection cuts, at the cost of a reduced event rate. Because the cross sections for the 
radiative top decay processes decrease whereas the itj rate increases when q top = —4/3, tt'y 
production is the largest contribution to the differential cross section over the entire Pr{l) 
range considered at both the Tevatron and the LHC (Figs. [5|b and f|b. The differential cross 
section of the t Why — > jjb'y "background" drops much faster with pri'y) than those 
of the two other contributions for large photon transverse momenta. When t — > Why — > 
jjbj selection cuts are imposed, qualitatively similar results are obtained, with the roles 
of the t — > Wb'y — > ivby and t — > Wb'y — > jjb'y curves exchanged and somewhat reduced 
contributions from tt'y production. 

The preceding results are summarized in Figs. ^ and [| Fig. |5|a shows the total SM 
(Qtop = +2/3) cross section for pp — > r yi ± vbbjj at the Tevatron, with all contributions from 
radiation in top production and radiative decays combined. The three curves correspond to 
the three sets of selection cuts. Fig. [|a shows the corresponding results for the LHC (with 
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pp initial state). We see in these figures, as we have already seen in Table |, that the ity 
selection cuts (solid curves) result in higher cross sections at both colliders than the cuts 
for radiative decays t — > Wiry. In Figs. |5|a and |6]a we also see that the photon distribution 
imposing t — > Wiry — > ivlry and t — > Why — > jjby selection cuts drops much faster than 
the Pt(t) distribution in the phase space region defined by the itj selection cuts. This is 
a consequence of the phase space restriction on the photon energy in radiative top quark 
decays. 

In Figs. |b and ||b we compare cross sections for q top = +2/3 (solid curves) and q top = 
—4/3 (dashed curves) for the Tevatron and LHC, respectively. In each figure the upper 
pair of curves corresponds to tt'y selection cuts and the lower pair includes events from 
both radiative decay regions combined. The figures show that at both machines, the effect 
of increasing \qt op \ is to increase the cross section in the tt'y cuts region and to decrease 
it in the radiative decay regions, for all values of photon p?. At the Tevatron, however, 
since the contributions from radiative top quark decays are small if tt'y selection cuts are 
imposed, changing the top quark charge has little influence on the shape of the photon p? 
distribution in this region (upper curves in Fig. |5|b). On the other hand, at small values of 
Pt(7), tt'y production may contribute significantly to the cross section in the regions defined 
by the t — ► Why — > ivlry and t — ► Wb'y — ► jjby selection cuts (see Table |TJ) . The shape 
of the photon transverse momentum distribution in these regions thus may be sensitive to 
q top (lower curves in Fig. ||b). The curves at the LHC (Fig. ||b) also exhibit some shape 
dependence, albeit less pronounced than at the Tevatron. 

We can look for sensitivity to q top in other distributions as well. In t — > Wb'y, emission 
of the photon from the 6-quark line leads to a collinear enhancement at small values of the 
separation between the photon and the 6-quark, AR(j, b) rad , which depends on q b = qt op —Qw- 
Therefore, in radiative top decays, the shape of the AR(j, b) ra d distribution is expected to be 
sensitive to the charge of the top quark. In contrast, the photon and the 6-quark originating 
from the t — > Wb decay in it —>■ r yW + bW~b events are essentially uncorrelated. In tt'y 
events, the photon and both 6-quarks are largely uncorrelated. Unfortunately, on an event 
by event basis, it is difficult to determine whether the b or the 6 quark originates from the 
radiative top decay. However, due to the collinear enhancement for small values of AR in 
radiative top decays, 

Ai?( 7 , b) min = min[A J R( 7 , 6), AR{ 1: b)\ (19) 

should coincide with the photon - 6-quark separation in radiative top decays for small values 
of AR. Alternatively, the AR(j, b s f) distribution, where b s j is the 6 or 6-quark with the 
smaller transverse momentum 

p T (b sf ) =mm[p T (b),p T (b)} (20) 

may be considered. Since a 6-quark which radiates a photon looses energy, the b s f in 
it — > , yW + bW~b frequently originates from the radiative top decay. It turns out that the 
AR(y, b s f) distribution more closely resembles the true AR(j, b) ra d distribution. Figure 
shows a comparison of the AR(j, b) ra d and AR(^y, b s f) distributions at the Tevatron in the 
SM case, imposing t — > Why — > tvby selection cuts. In each case, the individual distributions 
of the three contributing subprocesses are shown. Qualitatively similar results are obtained 
for the LHC case, and for t — > 1^67 — ► jjby selection cuts. Due to the collinear enhancement 
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at small opening angles between the 6-quark and the photon, the AR(j, b) ra d distribution 
rises very quickly with decreasing values of AR(j, b) ra d for the radiative top decay processes. 
For it^y production no collinear enhancement occurs for either the AR(j, b) or the Ai2(7,6) 
distribution. In Fig. 0a, we have chosen the AR(j, b) distribution for it^y production (dotted 
line). Very similar results are obtained if the AR(^y, b) distribution is chosen instead. Except 
for the steepness of the distribution at small AR(j, b s f) for t — > Wbq — > £^67, the shapes of 
the Ai?(7, b) ra d and AR(y, b s f) distributions are found to agree reasonably well. 

In Fig. [7| we have increased the minimum photon transverse momentum requirement 
of Eq. fli"3] ) to pr(j) > 15 GeV. As we have mentioned before, ity production contributes 
significantly to the cross section at the Tevatron in the region defined by the t — > Why — > 
Ivby selection cuts for the cuts specified in Sec. II. The ity contribution can be reduced by 
increasing the pt cut on the photon (see Fig. |3|), at the cost of a decreased number of events. 
Increasing the photon transverse momentum cut to Pr{l) > 15 GeV reduces the ity cross 
section for t — > Wby — > tvby selection cuts by about a factor 2, whereas the contributions 
involving radiative top decays decrease by only 30%. 

We compare the AR distributions for the two top quark charges in Figures [8] and |[ We 
display the normalized distribution (1/cr) (da/dAR(y, b s f)) for it — > •yivjjbb with q top = 2/3 
and q top = —4/3 at the Tevatron and LHC, respectively. Normalizing the distributions to 
unit area emphasizes the dependence of the shape of the AR(y, b s f) distribution to the top 
quark charge. For a charge —4/3 quark, the AR(j, b s f) distribution rises relatively faster in 
the region AR(j, b s f) < 1, but is flatter for 1 < Ai2(7, b s f) < 3, than if q top = 2/3. At small 
values of AR(y, b s f), emission from the 6-quark line dominates. Since the magnitude of <#> is 
the same for q top = 2/3 and q top = —4/3, the differential cross section is similar in both cases. 
However, the total cross section for q top = —4/3 within cuts is significantly smaller than for 
a SM top quark charge assignment (see Table |), due to destructive interference between the 
contributing Feynman diagrams. This leads to a larger normalized differential cross section 
for q top = —4/3 at small AR(j,b s f). In the region 1 < AR(j,b s f) < 3, the destructive 
interference is responsible for the flatness of the (1/cr) (da j 'dARi^f ,b s f)) distribution in the 
Qtop — —4/3 case. For Ai?(7,6 s j) > 3, the 6-quark with the smaller p T most of the time 
coincides with the bottom quark from the (regular) t —> Wb decay and, thus, carries little 
information on q top . The shape of the Ai?(7, b s f) distribution in the phase space region 
defined by the t — > Wb^y — > jjb r y selection cuts (Figs, and ||b) is found to be somewhat 
more sensitive to qt op than in the case when t — > Wb~f — > ivtry selection cuts are imposed 
(Figs. |8|a and ^a). At the Tevatron (Fig. |8|), we have again imposed a Pt(i) > 15 GeV cut. 
The shape of the AR(j, b s f) distribution for it — > r yW + bW~b at the LHC is in general more 
sensitive to the top quark charge than at the Tevatron. 



IV. QUANTITATIVE ESTIMATES 

In this section we determine quantitative bounds on 

2 , 
Aq top = qtop - - (21) 

which one expects to achieve with 20 fb _1 of data at the Tevatron in Run II, and 10 fb _1 at 
the LHC, analyzing the photon px and the AR(p{, b s f) distributions. We also briefly consider 
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the ratio of the cross section associated with radiation in top production to the radiative 
decay cross section which may be useful for small data sets. 

In our analysis of differential cross sections, we calculate 95% confidence level (C.L.) 
limits performing a x 2 test- The statistical significance is calculated by splitting the selected 
distribution into a number of bins, each with typically more than five events. In each bin 
the Poisson statistics are approximated by a Gaussian distribution. We use the photon 
transverse momentum distribution in the phase space region defined by the itj selection 
cuts of ([RJ), and the pr{l) and AR(j, b s f) distributions of the combined radiative top decay 
regions, defined by requiring that events satisfy QT7D or fllgp but not flTBp. We impose the 
cuts described in Sec. II, except for the AR(j, b s f) distribution, where we replace the cut 
on the photon of Eq. (jl3|) with Pt{i) > 15 GeV. We combine channels with electrons and 
muons (both charge combinations) in the final state. We also make use of the approximate 
symmetry of the cross sections under the transformation (§) and define the 6-quark charge 
to be qi, = q top — 1. Confidence levels calculated for q top = +4/3 will thus automatically 
apply for a charge —4/3 quark which decays into a 6-quark and a W~. 

Our calculation does not include QCD corrections which are expected to modify the cross 
section by 20 — 40%. In order to derive realistic limits, we therefore allow for a normalization 
uncertainty of AW = 30% of the SM cross section. The expression for \ 2 which is then used 
to compute confidence levels is given by [14j] 

* = Z®^H»u-i), (22) 

i=l J^i 

where rip is the number of bins, iVj is the number of events for a given Aq top , and N® is the 
number of events in the SM in the ith bin. / reflects the uncertainty in the normalization 
of the SM cross section within the allowed range, and is determined by minimizing \ 2 '- 

f (1 + AM)- 1 for f < (1 + AW)~* 
f=\f for (1 + AATp 1 < / < 1 + AAf (23) 

[ 1 + AW for / > 1 + AW 

with 

{no ~\ _1 n D aj-2 
£*? El- (24) 
i=l ) j=l iv i 

In Table we display the resulting 95% confidence level limits for Aq top expected from 
the Tevatron Run II for an integrated luminosity of 20 fb _1 , and from the LHC for 10 fb _1 . 
The last row shows the combined limits from the three individual distributions studied here. 
Possible backgrounds are assumed to be subtracted. At both the Tevatron and the LHC, the 
Pt(i) and AR(j, b s f) distributions in the combined radiative top decay regions yield rather 
similar limits for Aq top . As mentioned before, it'j production at the Tevatron is completely 
dominated by qq annihilation and photon radiation from the initial state quarks limits the 
sensitivity of the cross section to q top - Therefore, the bounds which can be achieved at the 
Tevatron from analyzing the photon pt distribution in the phase space region defined by 
the it'j selection cuts are significantly weaker than those which can be obtained from either 
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the Pt{i) or the AR(j, b s f) distribution in the combined radiative top decay regions.^] At 
the LHC, most of the it'y cross section originates from gluon fusion, and the limits from the 
photon pt distribution for it'j selection cuts are similar to those which can be achieved in 
the phase space regions where radiative top decays dominate. 

While a precise measurement of the top quark charge will not be possible in Run II of 
the Tevatron, even with an integrated luminosity of 20 fb _1 , we note that a charge —4/3 
top quark can be ruled out at « 95% CL. when the limits from the different distributions 
and phase space regions are combined. At the LHC, with 10 fb" 1 , a measurement of the 
top quark charge with an accuracy of about 10% appears to be feasible. As we pointed out 
in Sec. II, the pt cuts we use are suitable for running at low luminosity at the LHC, i.e. for 
C = 10 33 cm _2 s _1 . At high luminosity, £ = 10 34 cm~ 2 s~ 1 , the transverse momentum cuts 
must be increased due to pile-up effects caused by the large number of interactions per beam 
crossing. The , y£ ± i>bbjj cross section is quite sensitive to these cuts. For example, increasing 
the pt cut of the 6-quarks to 30 GeV and the photon, jet and missing p T cuts to 50 GeV 
reduces the ££7 and it — > r yW + bW~h cross sections by a factor m 10 and m 100, respectively. 
It is thus not clear whether running at high luminosity will improve the measurement of the 
top quark electric charge at the LHC. 

The bounds on Aq top are in general quite sensitive to AJV. If the normalization uncer- 
tainty can be reduced to AAf = 10%, the limits listed in Table [TI| can be improved by up 
to a factor 2. Further improvements may result from using more powerful statistical tools 
than the simple \ 2 test we performed. 

For integrated luminosities smaller than 20 fb _1 at the Tevatron, the limited number 
of events makes a detailed analysis of differential cross sections difficult. In this case, the 
ratio of the cross section associated with radiation in top production to the radiative decay 
cross section, 1Z, may be more useful. Many experimental and theoretical uncertainties, 
for example, those associated with the particle detection efficiencies, the uncertainty in 
the integrated luminosity, or the uncertainties from the choice of scale or of the parton 
distribution functions, are expected to cancel, at least partially, in 1Z. The cross section 
ratio is quite sensitive to q top . From the individual cross sections listed in Table |I|one finds 
1Z = 1.73 for q top = 2/3, and 71 = 2.88 for q top = —4/3 at the Tevatron. For an integrated 
luminosity of 10 fb -1 , the expected statistical uncertainty at the Tevatron is AlZ sta t = 0.65 
and a charge —4/3 top quark can thus be ruled out at ~ 90% C.L. using the cross section 
ratio, if one ignores possible systematic uncertainties. At the LHC, the value of 1Z is even 
more sensitive to q top than at the Tevatron, but enough events are expected to allow a more 
detailed analysis. 

V. SUMMARY AND CONCLUSIONS 

Like most of its fundamental quantum numbers, the electric charge of the top quark 
has not been measured so far. Alternative interpretations for the particle we believe is the 



3 In practice, hadronization effects we have omitted will come into play, and are likely to affect 
the AR distributions more than those for the photon pt- 
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charge 2/3 isospin partner of the 6-quark are thus not ruled out. For example, since the 
correlations of the 6-quarks and the VT-bosons in pp — > it — > W + W~bb are not determined, 
it is conceivable that the "t-quark" observed at the Tevatron is an exotic quark, Q4, with 
charge —4/3 which decays via Q4 — > W~b. This interpretation is consistent with current 
precision electroweak data [[|. In order to determine the charge of the top quark, one can 
either measure the charge of the 6-jet, or investigate photon radiation in it events. The 
latter method actually measures a combination of the EM coupling strength and the charge 
quantum number. Combining the results of the two methods will thus make it possible to 
determine both quantities. 

In this paper we have explored the possibility of measuring q top in pp — > r y£ ± b , bbjj, 
assuming that both 6-quarks are tagged. Our analysis makes use of both ity production and 
radiative top quark decays, and is carried out at the tree level. The EM coupling strength 
is assumed to take its SM value. Top quark and W boson decays are treated in the narrow 
width approximation. Contributions from radiative W decays can be suppressed by simple 
phase space cuts (see Eq. @). The remaining event sample can be separated into a ity 
sample, a it, t — ► Why — > yivb sample, and a it, t — > Why — ■> yjjb sample by imposing 
invariant mass and cluster transverse mass cuts (see Eqs. fll6|) - (|18D ). 

In the phase space region dominated by itj production, the cross section increases if 
the magnitude of q top increases. At the Tevatron, qq annihilation dominates and photon 
radiation off the incoming quarks constitutes an irreducible background which limits the 
sensitivity to q top - In contrast, at the LHC, gluon fusion dominates, and the ity cross 
section scales approximately with q\ . 

In the phase space regions which enhance it production with one of the top quarks 
decaying radiatively (t — > Wiry), interference effects between the Feynman diagrams where 
the photon is emitted from the top quark, the W boson and the 6-quark line may cause 
the cross section to either decrease or increase with Aq top = qtop — 2/3. For example, for 
qtop = —4/3, the yi^ubbjj cross section via radiative top quark decays is reduced by a factor 
of 3 to 4. 

The changes in the individual cross sections are reflected in the photon transverse mo- 
mentum distribution in the various phase space regions. We also found that the AR(y, b s f) 
distribution, where b s f is the b or 6-quark with the smaller p?, is sensitive to the charge of 
the top quark in the phase space regions dominated by it production where one of the top 
quarks decays radiatively. 

In order to determine how well one can hope to measure the top quark charge using 
photon radiation in top quark events in future Tevatron and LHC experiments, we have 
performed \ 2 tests of the Pt(i) and the AR(y, b s f) distributions in the various phase space 
regions of interest. At the Tevatron, with an integrated luminosity of 20 fb _1 , one will be 
able to exclude at « 95% C.L. the possibility that an exotic quark Q4 with charge —4/3 
and not the SM top quark has been found in Run I. For smaller integrated luminosities, the 
number of events expected is very small and it will be difficult to perform a quantitative 
analysis using differential cross sections. In this case, the ratio of the cross section associated 
with radiation in top production to the radiative decay cross section may be a useful tool. 
At the LHC with 10 fb" 1 obtained at C = 10 33 cm s , it should be possible to measure 
the electric charge of the top quark with an accuracy of about 10%. For comparison, at a 
linear collider with y/s = 500 GeV and JCdt = 200 fb _1 , one expects that q top can also be 
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measured with a precision of about 10% [15[]. Finally, at a 77 collider it is conceivable that 



the top quark charge can be determined with an accuracy of better than 1%, if the 77 — > it 
cross section can be measured with a precision of 2% [|16|]. 

The measurement of the electric charge of the top quark at the Tevatron using photon 
radiation in top quark events is severely limited by statistics. Even with an integrated 
luminosity of 20 fb _1 , qt op can only be determined with a precision of 30 — 100%. However, 
as we pointed out in Sec. I, information on the electric charge of the top quark may also be 
obtained from a measurement of the charge of the 6-jets, and the charge sign of the leptons 
in semileptonically tagged 6-quarks. This may significantly improve the precision which can 
be obtained. Detailed simulations will be necessary in order to determine how well the top 
quark charge can be measured using these methods. 
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TABLES 



TABLE I. Integrated cross sections for -yi vbbjj production at the Tevatron (pp collisions at 
= 2 TeV) and the LHC (pp collisions at yfs = 14 TeV) for q top = 2/3 and qt op = —4/3 imposing 

a) ity selection cuts, b) t — > Wby — > £u1ry selection cuts, and c) t — ► Why — ► jjby selection cuts. 

In addition, the cuts described in Sec. II are imposed. In each case, the contributions from it"f 

production, it production with t — ► Why — > iuby, and it production with t — ► IE67 — > jj'67 are 

shown. 



a) tt7 selection cuts 






«7 


t - 


-> IE&7 -> ^67 t 


— ► Why — > jj'67 








contribution 




contribution 


contribution 


total 


Tevatron, q top 


^3 


1.70 fb 




0.12 fb 


0.10 fb 


1.92 fb 


Tevatron, gt op 


1 

3 


2.37 fb 




0.03 fb 


0.02 fb 


2.42 fb 


LHC, qtop = 


^3 


33.3 fb 




5.78 fb 


4.86 fb 


44.0 fb 


LHC, qtop = 


4 
3 


111.4 fb 




1.21 fb 


0.97 fb 


113.6 fb 






b) t-> 


Why — ► £1/67 selection cuts 






Tevatron, q top 


^3 


0.36 fb 




0.26 fb 


0.13 fb 


0.75 fb 


Tevatron, q top 


4 

~~ 3 


0.52 fb 




0.07 fb 


0.03 fb 


0.62 fb 


LHC, qtop = 


^3 


0.68 fb 




3.67 fb 


2.26 fb 


6.61 fb 


LHC, qtop = 


4 
3 


2.25 fb 




0.87 fb 


0.59 fb 


3.71 fb 






c) t -» 


IE67 - 


> jjb'y selection cuts 






Tevatron, q top 


^3 


0.11 fb 




0.06 fb 


0.19 fb 


0.36 fb 


Tevatron, q top 


1 

3 


0.16 fb 




0.01 fb 


0.05 fb 


0.22 fb 


LHC, qtop = 


^3 


0.15 fb 




1.07 fb 


3.63 fb 


4.85 fb 


LHC, qtop = 


4 
3 


0.49 fb 




0.26 fb 


0.85 fb 


1.60 fb 



TABLE II. Limits achievable at 95% C.L. for Aqt op in pp — ► y^vbbjj at the Tevatron and 
the LHC. Bounds are shown for the photon transverse momentum distribution in the phase space 
region defined by the ity selection cuts, and the pt(i) and AR(-y, b s j) distributions of the combined 
radiative top decay regions. The last row displays the combined limits from the three distributions. 
The 6-quark charge is given by (ft = qtop — 1. The cuts imposed are described in Sec. II. For the 
AR(y, b s f) distribution, the pril) cut of Eq. (13) has been replaced by pr(l) > 15 GeV. 





distribution 


Tevatron, JCdt = 


= 20 fb^ 1 


LHC, JCdt = 10 fb" 1 


Pt(j), 
AR(y,b s 


prd)-, itj region 
comb, t — ► Wb-y regions 
f), comb, t — ► W 67 regions 


-1.16 < Aq top 
-0.32 < Aqt op 
-0.29 < Aqtop 


< 1.11 

< 1.20 

< 1.10 


-0.21 < Aqtop < 0.24 
-0.16 < Aqtop < 0.18 
-0.15 < Aqtop < 0.18 




combined 


-0.21 < Aq^ 


< 0.65 


-0.067 < Aq^p < 0.070 
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FIGURES 




FIG. 1. The photon transverse momentum distribution for pp — > 7^ vbbjj at the Tevatron, 
imposing ££7 selection cuts (see Eq. (JIB])), for a) (/top = 2/3 and b) qt op = —4/3. Shown are 
the itj (solid line), it, t — ► VF67 — > £1/67 (dashed line), and it, i — > VF67 — ► jj'67 (dotted line) 
contributions. The additional cuts imposed are described in Sec. II. 
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FIG. 2. The photon transverse momentum distribution for pp — > jt vbbjj at the LHC, impos- 
ing itj selection cuts (see Eq. (|Iq) ), for a) qtop = 2/3 and b) qtop = —4/3. Shown are the itj (solid 
line), it, t — ► Wiry — > tvb^ (dashed line), and it, t — ► Wb^f — > jj'67 (dotted line) contributions. The 
additional cuts imposed are described in Sec. II. 
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FIG. 3. The photon transverse momentum distribution for pp — > 7^ vbbjj at the Tevatron, 
imposing i — ► VF67 — > £^67 selection cuts (see text), for a) gt op = 2/3 and b) gt op = —4/3. Shown 
are the tit, t — > PT&7 — > £1/67 (solid line), tit, t — > ^67 — > jj'67 (dashed line) and itj (dotted line) 
contributions. The additional cuts imposed are described in Sec. II. 
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FIG. 4. The photon transverse momentum distribution for pp — ► j£ vbbjj at the LHC, im- 
posing t — ► ^67 — > £1/67 selection cuts (see text), for a) (/top = 2/3 and b) (ft p = —4/3. Shown 
are the it, t — > PT&7 — ► £1/67 (solid line), ft, t — > 1^67 — > jj'67 (dashed line) and itj (dotted line) 
contributions. The additional cuts imposed are described in Sec. II. 
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FIG. 5. The photon transverse momentum distribution for pp — ► j£ vbbjj at the Teva- 
tron. Part a) displays the individual contributions for qtop = 2/3, imposing itj (solid line), 
t — > — ► ^67 (dashed line) and t — > — > jj'67 selection cuts (dotted line). In part b), the 
photon distributions for the phase space region defined by the ttj selection cuts and the two 
radiative top decay regions combined are shown for q top = 2/3 and qtop = —4/3. The cuts imposed 
are described in the text. 
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FIG. 6. The photon transverse momentum distribution for pp — > j£ ubbjj at the LHC. Part a) 
displays the individual contributions for qtop = 2/3, imposing itj (solid line), t — > — > £1/67 

(dashed line) and i — > — > jj'67 selection cuts (dotted line). In part b), the photon 

distributions for the phase space region defined by the it r ) selection cuts and the two radiative top 
decay regions combined are shown for q top = 2/3 and qt op = —4/3. The cuts imposed are described 
in the text. 
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FIG. 7. Comparison of the SM AR(j, b) ra d and Ai?(7, b s f) distributions for pp — > ^Ivjjbb at 
the Tevatron, imposing t — > Wbj — > £^67 selection cuts (see text). Shown are the contributions 
from it production with t — ► ^67 — > £^67 (solid lines) and i — ► ^67 — > j.7'67 (dashed lines), and 
it r ) (dotted lines) production. The cuts described in Sec. II are imposed, except for the minimum 
photon transverse momentum which has been increased to 15 GeV. 
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FIG. 8. The normalized (l/c) (AR(j, b s f)) distribution for pp — ► ^lujjbb at the Tevatron, 
imposing a) t — > Wbj — > Ivfry selection cuts and b) t — > Wbj — > jjb^y selection cuts (see text). 
Shown are the results for a SM top quark charge assignment (solid lines) and for q top = —4/3 
(dashed lines). The cuts described in Sec. II are imposed, except for the minimum photon transverse 
momentum which has been increased to 15 GeV. 
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FIG. 9. The normalized (1/cx) (Ai?(7, b s f)) distribution for pp — ► ^ivjjbb at the LHC, imposing 
a) t — > Wirf — > Ivfry selection cuts and b) t — ► Wbj — > jj'67 selection cuts (see text). Shown are 
the results for a SM top quark charge assignment (solid lines) and for q top = — 4/3 (dashed lines). 
The additional cuts imposed are described in Sec. II. 
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